Bogoliubov theory of Feshbach molecules in the BEC-BCS crossover 
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We present the Bogoliubov theory for the Bose-Einstein condensation of Feshbach molecules in a 
balanced Fermi mixture. Because the Bogoliubov theory includes (Gaussian) fluctuations, we can 
in this manner accurately incorporate both the two-body and many-body aspects of the BEC-BCS 
crossover that occurs near a Feshbach resonance. We apply the theory in particular to the very 
broad Feshbach resonance in atomic ® Li at a magnetic field of Bo ~ 834 G and find good agreement 
with experiments in that case. The BEC-BCS crossover for more narrow Feshbach resonances is 
also discussed. 

PACS numbers: 03.75.-b,67.40.-w,39.25.-|-k 



I. INTRODUCTION 

Since the achievement of Bose-Einstein condensation 
(BEC) in an atomic Bose gas P, H 0, the field of ul- 
tracold gases has been one of the most active fields in 
physics. This is primarily due to the fact that these 
gases are very clean, accessible, and easy to manipu- 
late. For both theorists and experimentalists, ultracold 
atomic gases have therefore become an important play- 
ground for the study of a diverse set of quantum phe- 
nomena, that often in other guises also appear in very 
different condensed-matter systems. At present also de- 
generate atomic Fermi gases are available 3 and have 
been at the center of attention in the last two years. 
In particular, the superfluid Bardeen-Cooper-Schrieffer 
(BCS) phase has been created in two-component Fermi 
mixtures |E IE S B S To overcome the prob- 

lem of the small critical temperature, which depends 
in general exponentially on the interaction strength, all 
experiments have made use of a Feshbach resonance 
[llL Ha . Il3| . In this manner the interactions between the 
atoms can be strongly enhanced by an external magnetic 
bias field, giving ri se to the BEC-BCS crossover phe- 
nomenon |14L flSl Il6| . As a result of the atomic physics of 
the Feshbach resonance, the nature of the Cooper pairs 
in the BEC-BCS crossover is, however, not solely deter- 
mined by the interaction strength or scattering length, 
but in principle also depends on the width of the Fesh- 
bach resonance. In the limit of an infinitely broad res- 
onance, the properties of the gas can be derived from 
a single-channel theory that requires only the resonant 
scattering length as an experimental input. In general, 
however, a two-channel theory is needed. This is in par- 
ticular true for the description of the wave function of the 
Cooper pairs that plays an important role in the BEC- 
BCS crossover as we will show in this paper. 

In the crossover experiments with atomic ^Li, on which 
we focus mostly in this paper, the atomic gas is spin po- 
larized with respect to the electron spin of the atoms. 
In the event of a collision, the two atoms therefore ap- 
proach each other in a triplet state, which is called the 



open channel. Due to the hyperfine interaction, there is 
a nonzero coupling to a closed singlet channel, which has 
a different interatomic interaction potential that, most 
importantly for our purposes, contains a bound molecu- 
lar state. This bound molecular state we call the bare 
molecule from now on. Moreover, the triplet state has a 
magnetic moment, and the energy of the two incoming 
atoms can thus be tuned by a magnetic field due to the 
Zeeman effect. Near a Feshbach resonance, the threshold 
energy of the open channel is exactly tuned near a bound 
state of the closed channel and the nonzero coupling be- 
tween the two channels causes the scattering cross section 
between the atoms to be resonantly enhanced, in a simi- 
lar manner as in many collision experiments discovering 
new elementary particles in high-energy physics. 

In this paper we aim at presenting a theory of the BEC- 
BCS crossover that includes (Gaussian) fluctuations on 
top of the mean-field theory. Moreover, we include the 
above mentioned two-channel physics of the Feshbach 
resonance exactly and in this way improv e up on the var- 
ious theories [IIIilllllMIIlllllllillllllllll 
that have also been used to discuss different aspects of 
the BEC-BCS crossover observed experimentally. Some 
consequences of our theory were previously obtained in 
Ref. [23. Here we considerably expand upon that work 
by giving a detailed derivation of our Bogoliubov, or 
random-phase approximation (RPA) , theory and present 
new results in the BEC, unitarity, and BCS limits of the 
theory. The main physical idea behind the theory is that 
the entire BEC-BCS crossover can be seen as a Bose- 
Einstein condensation of pairs or dressed molecules, that 
in the extreme BEC limit are the bare molecules, and 
in the BCS limit the atomic Cooper pairs. The wave 
function of these dressed molecules is given by the linear 
superposition 



) = VZx,n(r) I closed) -I- \/l^^Xaa(r)|open),(l) 

where r is the interatomic distance. The parameter Z 
is the probability of the dressed molecule to be in the 
bare molecular state in the closed- channel of the Fesh- 
bach resonance. The spatial part of the bare molecular 
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wave function is denoted by Xm(r). The spin part of the 
bare molecule is equal to [closed), in agreement with the 
fact that it is a bound state in the closed channel. The 
open channel has the spin state jopen) and the spatial 
part of the wave function in this channel is Xaa(r)- The 
latter corresponds to the well-known wave function of the 
Cooper pairs in BCS theory. That this is a sensible way 
to describe the system at low temperatures, follows from 
the fact that it contains the correct physics in both ex- 
tremes. In the extreme BEC limit the dressed molecules 
are just bare molecules, i.e., Z ~ 1. In the BCS limit 
we have Z ~ and the dressed molecules correspond to 
the BCS Cooper pairs with the wave function Xaa(r)- In 
between these two limits the probability Z changes grad- 
ually and in this manner describes the smooth crossover 
from the BEC to the BCS limit. This important micro- 
scopic parameter has indeed been measured recently for 
^Li by Partridge et al. [l3| and was shortly thereafter 
calculated theoretically for that case in Ref . 28] . 

The paper is organized as follows. Because of its im- 
portance for the many-body theory that is of most in- 
terest to us, we first briefly discuss in Sec. II the atomic 
physics that is required to account for the resonant in- 
teractions between two atoms. In particular, we give a 
simple two-body derivation of the probability Z , that is 
then generalized in Sec. Ill when we discuss in detail the 
many-body Bogoliubov theory for the Bose-Einstein con- 
densation of the Feshbach molecules. In Sec. IV we apply 
the theory to the broad Feshbach resonance of ^Li at a 
magnetic field of Bo — 834 G. We here present results for 
the probability Z, the spectral function of the molecules, 
and the speed of sound in the gas throughout the BEC- 
BCS crossover. In Sec. V we briefly discuss the qualita- 
tive difference that would occur in these quantities if the 
Feshbach resonance in ^Li would be much narrower. In 
Sec. VI we consider the BEC, unitarity, and BCS hmits 
of the theory. We end our paper in Sec. VII with some 
discussion and conclusions. 



II. TWO-BODY PHYSICS AND BACKGROUND 
SCATTERING 

For several reasons it is most convenient to start with 
the discussion of the relevant two-body physics of the 
Feshbach resonance. First, it is a natural way to intro- 
duce the probability Z, that describes the exact bound 
state wave function of the two-body Hamiltonian. Sec- 
ond, it gives a quantitative description for the problem 
that works quite well in the BEC-limit, where there are 
no free atoms in the gas and the Bose-Einstein conden- 
sate of dressed molecules is well described by two-body 
physics. Third, it illustrates how to incorporate not only 
the resonant scattering, but also the background scatter- 
ing of the atoms, which turns out to be very important 
in the case of the extremely broad Feshbach resonance in 
^Li of interest to us. 

For two atoms, we can simplify the Feshbach problem 



by first splitting off the center-of-mass motion. In that 
manner we arrive at the time-independent two-channel 
Schrodinger equation 

(H,. Km\ \ _r,{m\ 

VKm Sb ) [Hn.) ) ^ ) ' 

where i/aa — —h^^^/m + Vbg(r) is the Hamiltonian of 
the relative motion in the presence of the background in- 
teraction Vhg{r), lipa) is the component of the two-body 
wave function in the open (triplet) channel, and IV'm) 
is the component of the two-body wave function in the 
closed (singlet) channel. The bare detuning from reso- 
nance is given hy Sb, and V^m is the coupling between 
the open and closed channels. We formally eliminate {ipa.) 
from Eq. |2Jl, by making use of the fact that the energy 
of the exact molecular bound state that we want to de- 
termine never lies in the spectrum of the atomic Hamil- 
tonian i?aa, and insert a complete set of states in the 
resulting expression. In the case of a large and positive 
background scattering length, there exists a bound state 
\iph) just below the threshold of the scattering continuum 
of the atomic Hamiltonian 77aa- Since this bound state 
may be relevant we include it explicitly as well. In this 
manner we ultimately find 

E -6b ^ (V'm|Vam-^r— Km|l/'m) 

= I (^(^m|K.„|V^[+^);^:^(V^l+Vam|^„.) 
+ (^m|Km|V'b)-^-^(^b|Km|^m), (3) 

where 2ek — h^'k'^ /m are the energies of the atomic scat- 
tering states IV'k^'')- We focus on the first term in the 
right-hand side of Eq. lO that is due to the continuum 
of scattering states, and leave the second term due to the 
bound state for later. 

We would like to use the pseudopotential approxima- 
tion for the matrix element ('!/'in|Kim|'0k^'')- Therefore, 
we rewrite the first term by using the definition of the 
two-body T-matrix KglV'k^^) = Tbg(2ek)|k), where |k) 
denotes the incoming plane-wave state associated with 
the scattering state I'i/'k^^)- We thus obtain 




For long wavelengths, the scattering properties of any 
interatomic potential become hard-core like, as the in- 
coming particles have too low energies to probe the 
structure of the interaction potential. The background 
scattering is then described by a single parameter, the 
background scattering length Obg. Moreover, the con- 
servation of probability then leads to the simply result 
(k'|rbg(2ek)|k) ~ 47rabg?i^/m(l -|-«fcabg). As a result, we 
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can introduce the coupling parameter g, that describes 
the coupHng between the two channels by means of 



(V''m|Km|'0] 



1 + ika 



bg 



(5) 



Hence the continuum contribution in the right-hand side 
of Eq. © is 



1 



(2^)3 l + (fcabg)2^-2ek' 

The zero-energy limit of this expression leads in Eq. ^ 
to a (finite) renormalization of the bare detuning (5s to 
the experimentally observable detuning S — A/i(i? — Bq) 
that is exactly zero on resonance. Here A/i denotes the 
magnetic-moment difference between the open and closed 
channels of the Feshbach resonance. For ^Li we have 
A/i ~ in terms of the Bohr magneton /is. Subtract- 
ing this limiting value we have that 




1 



1 
2^ 



(7) 



with r] = g'^m^^'^/ATrh^ |2gj. We have written fiT.^iE) 
to indicate that this term gives the self-energy of the 
molecules, which follows from the fact that the energy of 
the exact molecular bound state is found from 



E = S + hY.,n{E) 



(8) 



Now we wish to include also the bound state in the 
case of a large and positive background scattering length. 
For the moment we do not assume that the bound-state 
energy lies at —h'^/ma^^, but we initially use i?b = 

— h^/mal, where ah may be different from Obg due to 
effective-range effects. To determine the contribution of 
the bound state, we compare the normalizations of the 
wave functions of the bound and scattering states. We 
have in the pseudopotential approximation that 



4+)(r) = (l-r^ 



^b(r) = 



f 



o-r/a-t, 



(9) 
(10) 



which implies that for small separations between the 
atoms 

1 



'^b(r) = - 



27raggab 



(11) 



We can now use the results of Eqs. lfTT|l and ISJ to cal- 
culate the self-energy contribution of the bound state. 
Subtracting again the zero-energy limit to renormalize 
the bare detuning Sb we find now 



2h^/ma? 



bg 



1 



1 



Eh/E 
(12 



Using this result we can show that if £'b 
have 



ri^^{E) = Tj- 



-h^ /ma 



bg 



we 



1 



(13) 



Note that compared to Eq. (TJ the bound-state contri- 
bution has in the limit of a very large positive back- 
ground scattering length resulted in the simple replace- 
ment |abg| — > —flbg, which shows that the precise value 
of the energy of the bound state is all important in de- 
termining whether the combined effect of the scattering 
continuum and the bound state shifts the dressed molec- 
ular energy up or down. Using Eqs. ((SJ and (|12|l we 
can find the binding energy of the dressed molecule in 
the 2-body limit. We have applied this result to mea- 
surements of the binding energy of '^''K molecules 
and find excellent agreement with experiment as shown 
in Fig. n It is interesting that in this case the binding 
energy —Eh is substantially larger than /ma^^ and it 
is therefore a better approximation to use Eq. ((TJ for the 
molecular self-energy than Eq. H13|l . The use of the lat- 
ter self-energy would shift the molecular binding energies 
upward compared to the result with no background inter- 
actions, whereas agreement with the experimental results 
requires a downward shift. 
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FIG. 1: Binding energies for the dressed molecules as a func- 
tion of magnetic field. The solid line is given by the solution 
of Eqs. and H12|l . where we have used E\, = —8.9 MHz 
for the bound-state energy in the background potential |3l| . 
The dashed line is the solution in the absence of background 
interactions. The points are from the experiment of Moritz 
et al. 30]. The inset shows how the binding energy behaves 
as the energy of the bound state in the background potential 
is approached and an avoided crossing occurs. 

Now we can go back to the quantity Z that describes 
the closed-channel component of the dressed molecule 
as in Eq. QJ. While calculating the self-energy of 
the molecules, we have used a normalization such that 
("0111 IV'm) — 1- This implies that the normalization of 
the atomic part of the dressed molecular wave function 
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equals (1 — Z)/Z^ so we have the relation 

1 - Z 



Z 



(V'al^a: 



1 



[E - H^Y 



VamlV'n 



B=e„(5) 



dE 



(14) 



where em (5) is the energy of the exact dressed molecular 
state and the solution of the equation e^n = S + h'E,^[e^) . 
From this relation it follows that 



Z = 



1 



dE 



(15) 

E=e^{S) 

The probability Z of the dressed molecule to be in the 
closed-channel state was measured in the experiments of 
Jochim et al. 32] and Partridge et al. ^3 for ^Li. Us- 
ing the two-body results derived so far, we can get a 
very good agreement on the BEC-side of the resonance, 
as is shown in Fig. [3 At the two-body level, no stable 
molecules exist for positive detuning. Moreover, in that 
case Z becomes zero exactly at the resonance. Both fea- 
tures disappear at the many-body level as we discuss in 
detail in the next section. 
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A. Equation of state 

We start the development of the Bogoliubov theory 
by introducing the annihilation operator ^/'m(x) for the 
bare molecules, and the annihilation operators (x) and 
ipi (x) for the atoms in the two different hyperfine states 
I t) and I I), respectively. The Hamiltonian can then be 
expressed in second quantization as l35l| 



4m 



+ ,5B-2MW'm(x) (16) 



J2 fd^^u^)(- 

— t I V 



2m 



x)V't(x)^x(x) + Vj(x)^|(x)-i/'m(x)j 

dxV^|(x)V|(x)VT(x)^i(x) . (17) 



We use functional methods to calculate the physical 
quantities of interest, as it is simpler than the opera- 
tor formalism when going beyond mean-field theory. In 
the end, both the two-body physics, that we treated in 
the previous section, as well as the many-body BEC-BCS 
crossover physics are incorporated in our theory. This is 
achieved by first writing down the grand-canonical par- 
tition function of the gas as the functional integral 



Z^. 



, (18) 



where the Euclidean or imaginary-time action is 



dr< i/[V'm>m,C:V'<T] 



^ J dn '(p*{x,T)h-^i;^{x,T) 



+ I -!/;;„(x,r)?i— -0m(x,T) 



(19) 



FIG. 2: The two-body Z as a function of magnetic field for 
^Li, (a) without background interactions, (b) with a constant 
background scattering length Obg ~ — 1360ao equal to the 
background scattering length at resonance, and (c) with a 
magnetic field- dependent background scattering length as de- 
termined in Ref. 33]. The points are from the experiment of 
Partridge et al. 10]. 



III. BOGOLIUBOV THEORY OF MOLECULES 

In this section we give a detailed derivation of the Bo- 
goliubov or RPA theory of the BEC-BCS crossover. We 
discuss the problems that arise when trying to numeri- 
cally calculate the fluctuation corrections. We also con- 
sider the long-wavelength limit of the theory. 



We then rewrite the open-channel part of this action us- 
ing Green's function methods, where the Green's function 
G(x, t;x',t') can be seen as a 'matrix' with space and 
time 'indices'. We obtain 

S = JdTdX'tp'^{yi,T) (^^-^ - 2/i^'0m(x,T) 

- JdTdxdr'dx' [ri{^,r),ij^{x,T)]nG-' ?*(x'',t') 

(20) 

The coupling between the atoms and the molecules is 
contained in the atomic 2x2 Nambu-space matrix prop- 
agator G(x, r; x', r'), i.e., G = G[ip^, tp^n] depends on the 
molecular field V'm(x, r). This is a convenient rewrite, as 
it will allow us to integrate out the atoms using Gaus- 
sian integration. It also makes a clear connection with 
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the operator formaHsm, as the propagator is related to 
the expectation value of the imaginary-time Heisenberg 
operators in the following manner 



G(x,t;x',t') 



= - T 



^I(x',r'),VT(x',T')l )), (21) 



(x, t) 



where T is the time-ordering operator. The 2x2 matrix- 
structure allows us for including the coupling to the 
molecules and in particular immediately incorporates the 
anomalous expectation values of the atomic gas that arise 
when the molecules Bose-Einstein condense. 

It appears that we have neglected here the background 
interaction between the atoms. However, we reintro- 
duce the effects of the background scattering on the cou- 
pling g once we go back to momentum space, which, 
in agreement with Eq. ((SJ, results in the substitution 
g ff/(l + ikahg) where k is the magnitude of the rel- 
ative momentum between the two atoms. Diagrammati- 
cally this procedure is shown in Fig. |21 




ki+k2 



FIG. 3: Momentum-dependent atom-molecule coupling. The 
incoming molecular field is represented by the double lines and 
the single lines represent the outgoing atoms. The coupling 
g depends on the relative momentum k — (k2 — ki)/2 of the 
outgoing atoms in the case of a non-negligible background 
interaction. 

Our aim is to exactly integrate out the atoms and 
expand the result up to second order in the fluctua- 
tions of the molecular field around its expectation value 
(■0111 (x,t)), which we assume to be real from now on, 
and thus arrive at a Bogoliubov theory for the bare 
molecules. To make the connection with BCS-theory ex- 
plicit, we write (7/>m(x, r)) = A/g, perform the substi- 
tution ipni (x, r) ^ a /g + tpni (x, t) , and use the Dyson 
equation 



(22) 



with 




A 

2m 



X (5(x-x')(5(r-r'), 



U , 5^^-^") V(x-x')<5(r~r') 

n V3^m(x,T) y ^ ' > 



The reason for this procedure is that G~^ is now exactly 
the usual BCS propagator, which we know how to deal 
with exactly from BCS theory. The 'gap' A is in our 
case proportional to the expectation value of the molec- 
ular field, i.e., it represents the molecular BEC. Two 
atoms can be annihilated through the formation of a bare 
molecule in the condensate. Vice-versa, two atoms can 
also be created from the molecular condensate. Diagram- 
matically, the propagator can be expanded as is shown 
in Fig. II 



a) 
b) 



FIG. 4: Diagrammatic expansion of the normal and anoma- 
lous BCS propagators. The thick lines represent the full BCS 
propagator, the thin lines represent the bare propagator. In- 
teractions with the molecular Bose-Einstein condensate are 
denoted by the dashed lines, a) The normal BCS propagator, 
b) The anomalous propagator. 

The self-energy term E accounts for the molecular fluc- 
tuations, and as such depends on the molecular field, 
S = S[^/;j^, ■f/'m]- The action has become quadratic in the 
atomic fields, and so we can now formally integrate out 
the atoms by doing the Gaussian integral. We make use 
of the fact that G can be seen as a matrix with space-time 
'indices' and obtain 

Zgr = j d[V'm]'^['/'m] exp | - S'm [?/'m: V'ml/Ti 

+ Tr[ln(-G-MCV'm])] }, (24) 

where the trace is over Nambu space, and over the space 
and (imaginary) time coordinates. We next expand this 
exact result up to second-order in the molecular fluctua- 
tions by means of 

Tr[ln(-G-i)] = Tr[ln(-G-i(l - GaS))] 

1 



Tr[ln(-G, 



Tr[GaS + -GaSGaS 



(25) 



(23) 



Since the self-energy TiTi is linear in tpn-i this is indeed 
effectively an expansion in the molecular fluctuations. 
The first-order term, Tr[GaS], is required to be equal 
to (2/z — 5B){'4'n\ + V'm)/^ because we want to expand 
around the minimum of the effective action. Using the 
expressions for E and Ga, and going over to momentum 
space, the second-order term can be written in the zero- 
temperature limit as 
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Tr [GaSGaS] = ^ [V-'m(k, i^n), V'ml-k, -iw„)] 



i;2i(k,ia;„) S22(k,?w„) 



'(/'m(k, ZCJ„) 
V'ml-kj-j'^n) 



(26) 



For example, Sn contains two normal fcrmionic propagators, which can be taken from Eq. (|23|l . In momentum space, 
this ultimately leads to 



1 



13 J (27r)3 



i?icj„^+ e(k^) - ^ z;itj„'_ + e(k'_) - /i 



(27) 



Here we have also introduced the BCS dispersion ?iti>a(k) = [(e(k) — /i)^ + A^/(l + a^^k^)]^/^, the atomic dispersion 

e(k) = ?i^k^/2TO, and the notation k'j_ ~ k/2 ± k' and = n/2 ± n'. The coupling g is dressed by the background 
scattering through |g(k)p = g"^/ (l + k^a^g) as explained previously. The two propagators can be reduced using partial 
fraction decomposition, resulting in the expression 

ul{k')ulik'_) 



1 



dk' 



f3 J (27r)3 



+ 



?iWa(k^) + hLUii{k'_) — ihuj„ 

ulik'^)vlik',) 
?iaja(k^) — ?iWa(k'_) — ihujr, 

vlik'^)ulik'_) 
-7icja(k'^) + hLUii{k'_) — ihuj„ 

vl{k'^)vlik'_) 
-7iWa(k^) — ?ia;a(k'_) — ihtUn 

I 



iTlUJn'^ - ?iWa(k',_ 



;iWa(k+ 



?ia;a(k',_ 



1 

iflUJn' — ?iti'a(kL) 
1 

iflLOn' +flLOg_ik!_) 
1 

iflLOn' — 7iti'a(k'_) 
1 



(28) 



where Ua(k) and 
tors, obeying 



(k) are the usual BCS coherence fac- 



2,,^^ _ ?lWa(k) + e(k) - IJL 



<(k) - 



2nLj,{k) 

fiuja.{k) - e(k) + 

2nu^{k) ■ 



(29) 
(30) 



From the terms in the right-hand side of Eq. H28|) , we can 
read off the different processes that lead to a dressing of 
the bare molecules. The ladder diagram which this ex- 
pression corresponds to is shown in Fig. Eji. It describes 
two fermions, propagating with momentum k+ and k_. 
In the BCS-picture, there are indeed four possible pro- 



cesses that renormalize the molecules, due to the possibil- 
ity for the molecules to break up into two quasi-particles 
(proportional to ulu^, a quasi-hole and a quasi-particle 
(proportional to u^v^), or two quasi-holes (proportional 
to vlvl). 

Next, we take the zero-temperature limit by replacing 
the sum over internal Matsubara frequencies n' by an 
integral. As a result, only two processes survive: the cre- 
ation of a quasi-particle pair, and the creation of a quasi- 
hole pair. Finally, we perform the same renormalization 
of the detuning as explained in Sec. II. In this manner we 
are ensured that ^ = corresponds also in the many-body 
theory exactly to the location of the Feshbach resonance. 
To conclude, the self-energy contributions are now 



/dk 
(2^ 



dk' ul{k\)ul{\^_) vl{k'^)vl{)^_) 1 

(27r)3 '-"^^ ^' I ihujn - riu^{k'+) - hu^{k'^) ihujn + riu^{k'+) + hu^{k'_) 2e(k') | ' 



hT,i2{k,lLUn) = 2y — — ^ [^(k )| < Ua(k+)Wa(k+)Ua(k_)Wa(k_ 



(2^)3 " [ +' +' -' -'{hLO,{k'+)+hLU,{k'^)r + {hLUr,)^j ' 

hT,22(k,iujn) = hT,ii{k,-iu!n), 

7iS2i(k, iw„) = hT,i2{k,-iuJri)- (31) 



In our calculations, we also take into account that the background scattering length abg(-B) and the coupling 
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g{B) depend on the magnetic field, which turns out to 
be very important for the extremely broad Feshbach res- 
onance used in all experiments with ^Li 33]. 

Diagrammatically, Fig. |31 illustrates the way the BCS 
propagators are included into the theory. The ladder 
diagram in a) is the normal self-energy ?iEii. The BCS- 
propagators can be expanded as is shown in Fig.01 which 
gives rise to contributions such as the one in c). The di- 
agram in b) is the anomalous self-energy h'Ei2. A con- 
tribution using bare propagators is shown in d). This 
shows explicitly that there is only an anomalous contri- 
bution for nonzero A. 



to the modified BCS gap equation 




FIG. 5: Diagrammatic representation of the molecular self- 
energy diagrams. The bare fermionic propagators are dressed 
through interactions with the molecular condensate. Single 
lines denote fermions, double lines denote molecules. For the 
fermions, the thin lines represent the bare propagator, thick 
lines represent the BCS propagator, as in Fig.|l| Interactions 
with the condensate are denoted by the vertices containing a 
dashed line. The effect of the background scattering on the 
vertices is not drawn explicitly, a) The normal self-energy 
REii. b) The anomalous self-energy ?iEi2. c) An example of 
a contribution to the normal self-energy REn. d) An example 
of a contribution to the anomalous self-energy ?i,Ei2. 
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l.9(k)|^ 
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2nujJk) 



2e(k) 



(32) 



As expected, the same gap equation can also be ob- 
tained from the requirement that the linear terms {Sb — 
2^)(V-'m+^,*„)+Tr[G'a?iE] in the effective molecular action 
vanish. 

Integrating out the atoms has led to an effective Bo- 
goliubov action for the bare molecules. The interaction 
effects between the molecules and atoms are contained 
in the normal and anomalous self-energies fiHij. Be- 
cause we have expanded in the molecular fluctuations 
only, the BCS propagators for the atomic fields are in- 
cluded everywhere, and the BCS state is fully incorpo- 
rated in the physics. Beyond that, we have constructed 
an effective action that also describes the fluctuations or 
noncondensed molecules, which serves to draw not only 
conclusions about the condensate depletion but also on 
the nature of the dressed molecules in the condensate as 
we will see in a moment. 

The inverse propagator of the molecules takes the form 



-riG- 



K>l2i —hG 



-1 _1_ 

22 ^ 



?iS2 



, (33) 



where 



-hG^^{k,iLj) 
-;iG^2^(k,zcj) 



ihuj + f.{'k)/2 



2/i, 
2/i 



(34) 



are taken from Eq. (|20|) . The molecular propagator 
can be written as the expectation value of time-ordered 
molecular creation and annihilation operators, similarly 
as the atomic propagator in Eq. H21() . Using these ex- 
pressions we can obtain the grand-canonical partition 
function Zg^ and the thermodynamic potential fl = 
— ln(Zgi.)//5 in the Bogoliubov approximation by per- 
forming the Gaussian integral over the molecular fluc- 
tuations. Ultimately we obtain in this manner 



Zg, =exp \ Tr [in(-G-i)] - -Tr [H-G-')] 



P{S^2^i)V^ 



(35) 



In terms of the above self-energies, the bare molec- 
ular condensate density A/g is determined by the ex- 
act Hugenholtz-Pines relation 2/1 = 6 + ?il]ii(0,0) — 
?iSi2(O,0) ,^J, which in our case turns out to be equal 



where V is the volume of the gas. As the molecular and 
atomic contributions have now mixed, it is appropriate 
to determine the equation of state by taking the deriva- 
tive of the thermodynamic potential with respect to the 
chemical potential. Hence, the total atomic density be- 
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h(3V 
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-Tr 



G, 



?2 
9/1 



Tr [G„ 



(36) 



Here the third PauU matrix in Nambu space is denoted 
by CT3. Because we are working at zero temperature and 
in the thermodynamic hmit, we can rewrite the trace as 
an integral over momenta and frequencies. In detail this 
means that Y.\,,n /^PV = I dc^dk /(27r)'*. 

The four terms in the right-hand side of Eq. H36|) have 
the following physical interpretation. The first term is 
the usual BCS expression for the density of atoms in the 
BCS ground state. The second term gives the contribu- 
tion of the bare molecular Bosc-Einstein condensate to 
the total atomic density, and the third term gives the 
contribution of the bare molecules that are not Bose- 
Einstein condensed, i.e., it describes the depletion of the 
bare molecular Bose-Einstein condensate. The fourth 
and last term is most difficult to understand and can 
be best explained by reformulating the equation of state 
in terms of dressed molecules instead of bare molecules. 
If the gas contains a Bose-Einstein condensate of dressed 
molecules, there is both a bare molecular contribution 
and an atomic contribution to the total atomic density of 
the Bose-Einstein condensate. These two contributions, 
together with the contribution of the unpaired atoms, are 
contained in the first two terms in the right-hand side of 
Eq. (|36l) . These two terms correspond to the mean- field 
theory of the BEC-BCS crossover. The third and the 
fourth terms in the right-hand side of Eq. (|36|l repre- 
sent the effects of fluctuations. Together they give es- 
sentially the bare molecular contribution and the atomic 
contribution to the total atomic density of the dressed 
molecules that are not Bose-Einstein condensed. In the 
Bogoliubov theory, therefore, the gas consists of unpaired 
atoms, Bose-Einstein condensed dressed molecules, and 
dressed molecules that are not Bose-Einstein condensed. 



B. Numerical methods 

Though the analytical expressions for our theory are 
quite compact and their behavior can be studied in all 
limits, solving the gap equation in Eq. 132|l and the equa- 
tion of state in Eq. (|36|l simultaneously is quite involved. 
It is important to examine these equations carefully, in 
order to avoid numerical difficulties. It should be noted 
that the integrand of the self-energy not only depends on 
the magnitudes of the internal and external momenta, 
but also on the angle between them, so the integration 



over internal momenta is in fact a two-dimensional one. 
In combination with the integration over the external 
momenta and Matsubara frequencies and the iteration 
procedure required to find a self-consistent solution at 
all magnetic fields, this leads to the practical necessity 
of choosing a fast method of calculation where possible. 
Convergence issues turn also out to be a problem, espe- 
cially in the case of a broad resonance. 

Our general procedure is to solve the gap equation and 
the equation of state iteratively, to find the desired value 
of the total atomic density. The input for this procedure 
is the gap A, which is completely fixed during the itera- 
tions. Then we take a chemical potential, and solve the 
gap equation to find the detuning. Finally, the equation 
of state gives the atomic density that we can compare 
with the desired experimental density. We adjust the 
chemical potential appropriately, and start the next it- 
eration. We perform this procedure for different A to 
obtain a self-consistent solution for various values of the 
magnetic field. The reason why we keep A fixed rather 
than the detuning, is that this ensures that a physical 
solution of the gap equation can always be found. 

We have found that the procedure of tracing over the 
Matsubara frequencies is very subtle. The first possibility 
of doing this is to just replace the discrete sum by an in- 
tegral, due to the zero-temperature limit, and to perform 
a straightforward integration over the frequencies along 
the imaginary axis of the complex frequency plane. This 
is shown in Fig. by the solid curve. In this figure, the 
frequencies are shown as a complex energy z = ifiUn- Us- 
ing this method, we calculate the equal-time propagators 
^„ G(k, iujn) — ^(k) -1-1/2, and therefore a convergence 
factor is needed to obtain the correct occupation number 
iV(k). This is equivalent to closing the contour on the 
left with a large semicircle. 



Im[z] 










[ 



Re[z] 



FIG. 6: Tracing the Matsubara frequencies using various 
methods of contour integration. 

For the third and fourth terms in the right-hand side 
of Eq. H3t)|l this method has a significant disadvantage. 
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In particular, for the third term and for high external 
momenta, the integration over the imaginary axis, and 
the closing of the contour add to zero, and this means 
the calculation would be very sensitive to the precision 
with which these contributions can be calculated. 

An alternative is to perform a contour deformation, 
such that the path of integration lies infinitesimally below 
and above the real axis, as is shown by the dashed lines 
in Fig. IHl This is usually a desirable method, as there is 
no need to correct for the value of the equal-time prop- 
agators. The use of this method, however, necessitates 
knowledge of the analytic structure of the retarded prop- 
agator to accurately perform the numerical integration, 
which is not feasible to obtain from our complicated ex- 
pressions for the self-energy. Therefore, we prefer to use 
another method by integrating over some intermediate 
path, and to work with complex energies along the inte- 
gration. Any parabolic shape would seem very suitable, 
as it also takes away the need for closing the contour on 
the left side. Unfortunately, these curves have a strong 
positive and negative contribution over a wide range of 
energies, that would again have to cancel. 

For the third term in the right-hand side of Eq. (|36|l . 
the best approach, shown by the dash-dotted line in 
Fig. El is to go along the imaginary axis for a finite in- 
terval cr, and then to follow a special path given by 



zip) =e(k)/2 - 2^ 



- (p+l)Ve(k)/2-2^-icr 



prj 



(37) 



where p parametrizes the path and e(k) is the external 
momentum. This expression is chosen such that the two- 
body self-energy would give exactly zero along the path. 
Closing the contour still gives a correction of 



-Arg 

TT 



r/+ v/e(k)/2-2^-i(T 



(38) 



where Arg[z] gives the argument of the complex number 
z. This term becomes small for small a and large e(k). 
For the fourth term in the right-hand side of Eq. H36(l . 
such a procedure is not possible, and we resort to a 
parabolic path around the negative real axis. In this 
case, the convergence is better, as the integrand behaves 
as l/|z|2. 



C. Long- wavelength limit 

We finally study the behavior of the poles of the bare 
molecular propagator. The molecular propagator at zero 
momentum describes the closed-channel part of the Bose- 
Einstein condensate of dressed molecules and therefore 
has a pole at zero frequency. The strength of this pole is 
of great interest to us, as we will see in a moment. 

The behavior of the propagator for small k and lu can 



be studied by expanding around this point. We have 



?iEn(k, z) = nS^^) + + 



1. 



^(0,2)^2 



?iEi2(k,z) =nsf/^ + -Eir'z^+ 
Slf^e(k) + ..., 
where the following definitions have been used 
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(39) 
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(40) 



k=O,z=0 



k=O,z=0 



-^?ii;i2(k, z) 



k=O,z=0 



We can find analytic expressions for all of these coeffi- 
cients using Eq. H31|l . 

Now we introduce the spectral function, by Wick- 
rotating the molecular propagator to real energies, i.e., 
Gm.ii(k, ztj) Gm,ii(k, 0;+). The imaginary part of 
this expression gives the spectral function of the bare 
molecules. This spectral function gives us the wave- 
function overlap between a bare molecule and the two- 
atom elementary excitation of our many-body system. 
For comparison, in the limit of two-body physics, the 
spectral function can be defined in a similar way. It is 
given by 

Pm(k,t^) 



.-ilm 

TT 



huj+ - e(k) /2-5- hT.mifu^+ - e(k) /2) 



(41) 



with the self-energy given in Eq. (jT)). For negative de- 
tuning, this spectral function contains a delta peak at 
negative frequencies, corresponding to the dressed molec- 
ular state, and a tail for positive frequencies associated 
with the two-body scattering states. As the delta peak 
describes the overlap between the dressed and molecular 
state, its strength is exactly equal to Z. 

In the many-body case, the situation is similar. At zero 
momentum, the delta peak lies at zero energy, because 
it is the probability of creating a dressed molecule in the 
Bose-Einstein condensate. Its strength is equal to the 
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sought-after factor Z. We derive analytically that 
1 



/3m(k, W) 



Im [Gn(k,u;+)] 

ty.(0:0)|i/2^ 

12 I \ iij-.uj^hck) 




E {huj + Tick) , (42) 
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with = (l-?iEi"^ 
speed of sound is given by 
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(43) 



This shows that our theory has reduced to a renormal- 
ized Bogoliubov-theory for the bare molecules. The delta 
peak at positive energy corresponds to a particle-like ex- 
citation, and its strength is equal to the Bogoliubov am- 
plitude Z|uni(k)p. The delta peak at negative-energy 
corresponds to a hole-like excitation and has strength 
Z|wi„(k)p. Both of these are renormalized by an overall 
factor Z, but from their separate expressions this is not 
immediately clear. However, as mentioned before, we ex- 
pect the combined effect of the two delta peaks at zero 
momentum to result in a single delta peak in the spec- 
tral function with a strength Z . From the expression in 
Eq. H42() this is found to be the case, and the factor Z is 
found to be equal to 



Z = 



1 



v(o,i) 
^11 



(1 



.(0,1) 

^11 



(0,2) 
2 



.(0,2) 
^11 



(44) 



IV. BROAD RESONANCES 

In the case of a broad resonance, such as the one at 
834 G for ^Li, both the molecular density and Z become 
very small towards the unitarity limit (5 = 0. We will 
discuss the physics of the system using two figures that 
summarize the main results of our calculations. 

In Fig. [71 the spectral function for the bare molecules 
is given for the BEC-side and the BCS-side of the reso- 
nance. The top shows that the BEC-limit shows a struc- 
ture that resembles the two-body limit. The delta peak 
is the bare molecular condensate that by definition exists 
at zero frequency. The continuum of atomic scattering 
states causes a long tail on the right side of the figure. 
The distance between these two features is well approx- 
imated by the two-body bound state energy €^{5). The 
relevant energy scale in this limit is much bigger than the 
Fermi energy, due to the extreme broadness of the reso- 
nance, which is the reason that two-body physics works 
so well. 
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FIG. 7: Spectral functions of the bare Li2 molecules with 
zero momentum a) in the BEG limit at 571 Gauss and b) in 
the BCS limit at 892 Gauss of the crossover occurring near 
the broad Feshbach resonance of ®Li at 834 Gauss. The Fermi 
energy of the gas is 380 uK. 



In the bottom half of the same figure we present the 
spectral density in the BCS-limit that shows completely 
different physics. Again, the bare molecular condensate 
is at zero frequency, which is now also the position of 
twice the Fermi surface, and the two-quasiparticle con- 
tinuum shows up as particle excitations on the right, and 
hole excitations on the left. There are no fermionic ex- 
citations within 2A of twice the chemical potential, as is 
expected from BCS theory. Left and right we see also 
the peaks that are expected from RPA theory and physi- 
cally correspond to the long-lived modes associated with 
the fluctuations in the magnitude of the BCS gap A. 
A slight curving at these ends anticipates these peaks, 
which is due to the momentum dependence of the cou- 
pling constant g. 

In Fig.|Hlwe then show with the solid line the closed- 
channel fraction Z as a function of magnetic field. This 
is an improved version of the figure in Ref. [2^, which 
contained a numerical error in the determination of Z . A 
comparison is made to the experimental data of Partridge 
et al. |l3- For reference, the bare molecular condensate 
density fraction is also plotted. These quantities are not 
identical and their ratio diverges as 5 in the BCS-limit, 
as we will show later on. The inverse of this ratio is equal 
to the condensate fraction of dressed molecules, and it is 
plotted in the inset. Fluctuations become very important 
on the positive side of the Feshbach resonance, and their 
contribution to the total density becomes of the same 
order as the condensate density. 

In Fig. 1^1 we see the behavior of the chemical poten- 
tial as a function of the magnetic field. For low mag- 
netic fields, the chemical potential is well approximated 
by half the binding energy of the two-body bound pair 
tra{5)- This corresponds well to our picture of having 
a pure Bose-Einstein condensate of dressed molecules in 
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FIG. 8: The solid curve shows the probabihty Z for ^Li. The 
data points are from the experiment of Partridge et al. [Tof . 
The dashed curve shows for reference the fraction 2Zn„ic/n 
as a function of magnetic field. The Fermi energy h^kp/2m of 
the gas is 380 nK. In the inset the solid line shows the Bose- 
Einstein condensate fraction of dressed molecules 2nmc/n and 
the dashed line the contribution of the fluctuations to the total 
atomic density. The vertical lines indicate the magnetic flelds 
where ftpiarosl = 1- 
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FIG. 9: Chemical potential as a function of magnetic held for 
"Li and for the same parameters as in Fig.|S| At unitarity, we 
find /i ~ 0.51eF. The dashed line is given by eni((5)/2, where 
Em (5) is the energy of the dressed molecular state in two-body 
physics. 

this limit. At high magnetic field, the chemical poten- 
tial saturates and approaches the Fermi energy, as is ex- 
pected for a fcrmionic gas. At zero detuning, we have 
that /i ~ O.Slep, which is commonly written as {\ + [3)e-p. 
The result of /3 ~ —0.49 can then be compared to the 
Monte-Carlo resuh of /3 ~ -0.58 37, 38]. 

The speed of sound, given in Fig. 1101 has an inter- 
esting behavior, due to the interplay of the resonance, 
the background scattering, and the dependence on the 
magnetic field of the different parameters. Apart from 
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FIG. 10: Speed of sound in the case of "Li and for the same 
parameters as in Fig. |H| At unitarity, we find c ~ 0.43wf. 
The dashed line is the BEC-limit derived from Eq. 1491 1 using 
Eqs. H46fl . I|47fl . and H5()|l . The dash-dotted line is the speed 
of sound wf / \/3 of the Anderson-Bogoliubov mode. 

the crossover near resonance, we see around a mag- 
netic field of 550 G, another crossover taking place 
from the resonant-scattering dominated regime to the 
background-scattering dominated regime. In the extreme 
BEC-limit, the sound velocity goes to zero, as the gas 
approaches the limit of an ideal, noninteracting gas of 
molecules. On the right it approaches the BCS-like limit 
of the Anderson-Bogoliubov mode. In the unitarity limit, 
we find c ~ 0.43uf. This sound velocity can be com- 
pared to previous work js^ where a unitarity value of 
c ~ 0.37i;f was found. One of the important reasons for 
this difference is the inclusion of background scattering 
effects in our approach. 

V. NARROW RESONANCES 

We now consider a much narrower resonance, in the 
absence of background scattering. From an experimental 
viewpoint, this resonance with if' ~ e-p, where ep is the 
Fermi energy, is an extremely narrow one, but from a 
theoretical point of view it is intermediate. 

Fig.^Jshows the closed-channel component Z, and the 
bare molecular condensate density in this case. Again, 
the difference between Z and the bare condensate frac- 
tion is clearly visible. What is interesting here, is that 
the crossover takes place at a somewhat higher detun- 
ing. At zero detuning we still have a large fraction of 
Bose-Einstein condensed dressed molecules, as they are 
stabilized by the Fermi-sea, which blocks the dissociation 
of the condensed molecules. This shift is also observable 
in Fig. ^1 where the chemical potential is plotted as a 
function of detuning. The chemical potential follows to a 
large extent the energy level of the bare molecule through 
the resonance, and only distinctly saturates at positive 
detuning. At resonance, we now find a lower value of 
/i ~ O.lTep. This shows that universality is lost for nar- 
rower resonances. 

The speed of sound, shown in Fig. 1131 has a fairly 
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FIG. 11: The solid curve shows the probabiUty Z for a narrow 
resonance in the absence of background interactions. The 
dashed curve shows the fraction 2Znmc/w as a function of 
magnetic field. In the inset the solid line shows the Bose- 
Einstein condensate fraction of dressed molecules 2nmc/n and 
the dashed line the contribution of the fluctuations to the 
total atomic density. The vertical line shows the position of 
the resonance. 
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FIG. 13: Speed of sound in the case of a narrow resonance 
[rj^ = ep). The dashed curve shows the BEC-limit for a nar- 
row resonance given in Eq. 15411 . The dash-dotted line is the 
speed of sound vf/V^i of the Anderson-Bogoliubov mode. 
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FIG. 12: Chemical potential as function of detuning for a 
narrow resonance {rf = ep). At unitarity, we find /i ~ 0.17eF. 
The dashed line is given by em(<5) /2, where em(i5) is the energy 
of the dressed molecular state in two-body physics. 



straightforward crossover behavior, mainly due to the 
absence of background scattering. Again we retrieve the 
Anderson-BogoUubov mode on the right, and a vanishing 
speed of sound on the left. 



FIG. 14: Logarithmic plot of the coefficients in the Bogoliubov 
Hamiltonian. The coefficients are defined in Eqs. 1401 . The 
term that dominates on the BEC-side is a first-order deriva- 
tive in energy (E^j'^'), on the BCS-side it is a second-order 
derivative in energy (E^j''^'). 
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VI. VARIOUS LIMITS OF THE THEORY 

In this section we study the three notable limits of the 
crossover: the BEC-limit, the BCS-limit and the unitar- 
ity limit in more detail. Focusing on the microscopic 
physics and the long- wavelength limit, the crossover can 
be studied by looking at the behavior of the self-energies 
for small energies. The various coefficients of interest are 



In the BEC-limit, we expect to retrieve the results of 
two-body physics. The reason for this is that, as a func- 
tion of the detuning, the chemical potential quickly be- 
comes very negative, eliminating the effects of the Fermi 
sea. The diagonal terms in the many-body molecular 
self-energy reduce to the two-body self-energies, and the 
off-diagonal self-energies become less important. We de- 
fine the BEC-limit by |A| ^ To leading order we 
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derive from Eq. (|31|l 



?iEii(2;,k) 



+ e(k)/2 - 2Ai 



l + ^^^^V-^ + <k)/2-2/i 



(45) 



in agreement with Eq. ((7|). Furthermore, we wish to de- 
rive the expression for Z given in Eq. (|15|l . This follows 
directly from our general Eq. H44|) for Z as S12 is small 
compared to the diagonal entry En. Thus 



Z 



1 - 



= 1 



2V^(1 + ^^^^^y=2]I) 



(46) 



in agreement with Eq. p5|) . In the BEC limit the fluctua- 
tion terms in the equation of state can be neglected. This 
then suggests that we have only a Bose-Einstein conden- 
sate of dressed molecules, and nothing else. Our bare 
molecules are contained in the term 2A^/g^, the density 
of bare atoms is given by the term Tr[G'a(T3]. Indeed, 
defining the density of Bose-Einstein condensed dressed 
molecules as rime = f^/2, where n is the total density, 
these terms and the expression for Z reduce to 



Tr[Gaa3] = 2(1 - Z)n„ 



(47) 
(48) 



This is illustrative for the fact that the dressed molecules 
can be split in a molecular component with contribution 
and an atomic component with contribution \ — Z . 
The correct BEC limit for the speed of sound is slightly 
more involved, as it is not a result of two-body physics. 
It explicitly depends on the off-diagonal terms in the self- 
energy as we show next. We can make a straightforward 
simplification by again using A as our small parameter. 
Furthermore we use the fact that sj"'^'' = — 2I]^^'*''' in 
this limit. The general expression for the speed of sound 
is then 



2m ( 1 - ^' 



1/2 



(49) 



Depending on the background scattering length, the 
width of the resonance, and the density, the system can 
enter different regimes. Expanding in the parameter 
A/ 1^1, the off-diagonal term reduces to 
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6OV2 f -l^y + 80 f -^V^' + 20x/2-^ + ^ , 



(50) 



where i?bg = ^ /™ctbg- '^^^ S^-P is fundamentally a 
many-body property, but it can be related to the two- 
body Z and the density through Eq. (|47ll . One possible 



regime is ^ E'bg, where the background scattering 
is small, while being sufficiently close to resonance such 
that I (5 1 ^ 77^, yet sufficiently far to be in the BEC-limit. 
In that regime the resonant scattering determines the 
speed of sound, and wc find the simple expression 



/m. 



(51) 



where arcs = Ti-q/S^/rn ~ fir] / /S.^{B ~ Bq)^^ is the reso- 
nant scattering length, and 7ii„c is the Bose-Einstein con- 
densate density of dressed molecules. This can be cast in 
the more familiar result 



c = ?iV47rai„ninc/ "in 



(52) 



using the relevant parameters for molecules = 2ares 
|40l | and m„i — 2m. This is precisely the Bogoliubov 
speed of sound j41J. Note that the Bose-Einstein con- 
densate density of dressed molecules shows up in the ex- 
pression for the speed of sound, further strengthening the 
case that the dressed molecules are the relevant physical 
entities in this system. In the case of ^Li, at densities that 
are typical for experiment, this regime is however not well 
separated from other regimes, due to the extremely large 
background scattering. 

Before we move on to the background-scattering dom- 
inated regime, we compare our results to single-channel 
theory. In this case there is only a single parameter, the 
total atomic scattering length a = a^cs + Obg- Assuming 
the validity of Eq. 1)52(1 . the molecular scattering length 
can be derived from the speed of sound, and expressed 
in terms of the total atomic scattering length a. This is 
shown in Fig. 1151 The expected one-loop single-channel 
result am = 2a is given by the dashed line and can ba- 
sically never be considered as a good approximation to 
the one-loop result obtained from our two-channel theory. 
We expect the same to be true of the exact single-channel 
result. 




500 



550 



600 650 700 
Magnetic field (G) 



750 



800 



FIG. 15: Ratio of the molecular and total atomic scattering 
length as a function of magnetic field. 

When the background scattering gives a significant 
contribution, and ~ -E-bg, which for ^Li at the densi- 
ties of interest happens roughly around a magnetic field 
of 650 G, it has a strong effect on the speed of sound. 
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Taking it into account using Eq. (|50|l . and relating the 
speed of sound to the two-body Z as described above, we 
get the result that is drawn by the dashed line in Fig. 1101 
We can conclude that the speed of sound is not so well ap- 
proximated with two-body physics as the closed-channel 
component Z . This is maybe to be expected as the speed 
of sound is really a many-body property of the gas. 
In the extreme BEC-limit, the speed of sound vanishes 



3/2 



2[-em(<5)]= 



1/2 



(53) 



small A. We can use Eq. (|43|) to find the speed of sound 
in this limit and we obtain 



73' 



(56) 



where vp is the Fermi velocity fikY/m. This coincides 
with the speed of sound of the Anderson-Bogoliubov 
mode 0, 01 . A more illustrative way of deriving 
this result is to determine the eigenvalues of the prop- 
agator in the long- wavelength limit. The leading terms 
are the ones mentioned above, leading to 



where ein(^) is the energy of the exact bound molecu- 
lar state taken from Eq. ||SJ|. In this way the theory is 
completely reduced to two-body physics, where there in- 
deed exists no sound. Of course, the vanishing of the 
speed of sound in the extreme BEG limit is a result of 
the fact that we neglected the background scattering be- 
tween molecules. If this is effectively repulsive, it will 
ultimately stabilize the molecular condensate and result 
in a nonzero speed of sound. Finally, in theory we can 
also be far away from resonance, and in a regime where 
the effect of the atomic background scattering is small. 
In this case the speed of sound reduces to 



1/2 



(54) 



^8(-(5)3/2 

This solution is shown by the dashed line in Fig. El 

B. BCS-limit 

For large positive detuning, the system strongly resem- 
bles that of a normal BCS superfluid. The density of bare 
molecules becomes exponentially small, and the proper- 
ties of the system depend mainly on the gap parameter 
A, which is also decreasing exponentially. Moreover, the 
chemical potential tends to the Fermi energy. 

Again, we can study the system by looking at the ac- 
tion for long wavelengths. We find explicitly, 



-V 



2V2 

TT 



1 



2eF/-Bbg 



log(A/eF), 



2eF/-Ebg' 
2%/2 ri l-2eF/£;bg 



hi: 



(0,2) 
12 

(1,0) 
11 

(1,0) 
12 



TT 7?F4(l + 2eF/^bg 
-2V2e^77/37rA2 
\/2e^?7/37rA2 



■log(A/eF), 



(55) 



3/2 



r;/97rA2 



where ep is the Fermi energy. We see that the coefficients 



^(0,2) 
^11 



.(0,2) „(1,0) 



^12 



^11 



and 2^2 *^^ diverge most strongly for 



367rA2 



(57) 



6{hujf + SekEF -3(;iw)2 - 4ekeF 



-3(?ia 



4ekeF 6{huj)- + Se^ep 



-011 



This action can be diagonalized by introducing amplitude 
and phase fluctuations. Keeping in mind that ip^n are 
the fluctuations of the molecular field, we define ipm = 
p/g + i^d/g and V'm — P/d ~ to obtain 



187r.g2A2 



k,n 



9{hLjy + 12e(k)ep 
~^2^[ a5 +3|0(k,*w„)l 

k,n ^ 

: {3{hujf + 4e(k)eF) , 



P 

A9 



(58) 



in agreement with Ref. [43 • Here N{0) is the density 
of states of a single spin state at the Fermi surface. We 
expect that the non-leading terms in A lead to a mass 
term of order A^ for the amplitude fluctuations. The 
result in Eq. H56() can be immediately derived from the 
phase fluctuation part. 

The closed-channel component in the BCS-limit is 



Z = 



^(1 - 2ep/£;bg) 



3/2 



A2log(ep/A). 



(59) 



It is important to notice that in this BCS limit, Z and the 
density contribution of the Bose-Einstein condensate of 
bare molecules 2A'^/g'^ have an entirely different behavior 
as a function of detuning. From the gap equation it can 
be derived that the gap vanishes exponentially. 



A c>c ep exp < — 



^(1 + 2ep/£;bg) 



(60) 



thus the ratio between Z and A^ diverges linearly with 
the detuning. This shows that these are different entities 
with different physical meanings. 
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C. Unitarity limit 

For 6 — we can make a rough estimate of Z, that 
gives a simple physical reason for why Z is small for 
broad Feshbach resonances. As can be seen in Fig. [7| 
the most notable result of the many-body theory is the 
gap around the Fermi surface. In contrast to two-body 
physics, where Z goes to zero at resonance, the many- 
body physics stabilizes the pairs, much like Cooper-pairs. 
For broad resonances, the spectral function has a long 
tail at high energies, which is not affected by many-body 
physics because it lies far above the Fermi energy. Using 
the fact that the spectral function is normalized to one, 
we can thus approximate the spectral weight in the delta 
function by integrating the two-body spectral function 
over an energy interval of twice the gap. 

TT Jq Z"^ + rj''Z TTT] 

From our calculations, the gap at unitarity is equal to 
0.46eF, which is quite close to the Monte-Carlo result of 
A ~ 0.50eF (Hi- Substituting this in Eq. ^ results 
in an approximate value for Z of 3.5 • 10^^, which is 
only about 20% lower than the result of our Bogoliubov 
theory. 



VII. CONCLUSIONS 

The BEC-BCS crossover near a Feshbach resonance is 
very rich in physics. We have presented a theory that rig- 
orously incorporates the two-body Feshbach physics and, 
due to the inclusion of fluctuations, allows for an accu- 
rate study of the full crossover problem, both for broad 
and narrow resonances. We calculated the probability Z, 
that is crucial in understanding the microscopic physics 
of the crossover, and studied the molecular spectral func- 
tion and the spectrum of collective modes throughout the 
BEC-BCS crossover. Where possible, we have compared 
our theory with experiment and have found a satisfying 
agreement, even in the unitarity limit where the interac- 
tions are very strong. 
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